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Abstract. It has recently been suggested that ultra- high energy cosmic rays could 
have an extragalactic origin down to the "second knee" at ~ 4 x 10"'^^ eV. In this case 
the "ankle" or "dip" at ~ 5 x 10^* eV would be due to pair production of extragalactic 
protons on the cosmic microwave background which requires an injection spectrum of 
about E~^'^. It has been pointed out that for injection of a mixed composition of 
nuclei a harder injection spectrum ^ E~'^ '^ is required to fit the spectra at the highest 
energies and a galactic component is required in this case to fit the spectrum below 
the ankle, unless the proton fraction is larger than 85%. Here we perform numerical 
simulations and find that for sufficiently magnetized sources, observed spectra above 
10^^ eV approach again the case of pure proton injection due to increased path-lengths 
and more efficient photo-disintegration of nuclei around the sources. This decreases 
secondary fluxes at a given energy and thus requires injection spectra ~ E~^'^ , as 
steep as for pure proton injection. In addition, the ankle may again be sufficiently 
dominated by protons to be interpreted as a pair production dip. 
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1. Introduction 

A major unresolved aspect of ultra-high energy cosmic ray (UHECR) physics ^ is their 
composition and above which energy the flux is dominated by extragalactic sources. 
Above ~ 10^^ eV the chemical composition is basically unknown |2]. Around 10^^ eV the 
situation is particularly inconclusive as HiRes data suggest a light (proton dominated) 
composition, whereas other experiments indicate a heavy composition j3] . 

As a consequence, there are currently two different scenarios: The "standard" one, 
where a transition from a steeper, galactic heavy component to a flatter, extragalactic 
component dominated by protons takes place at the ankle at ~ 5 x 10^^ eV, see, e.g., 
Ref. ini, and a more recent one suggesting that this transition actually takes place at 
lower energies, namely around the "second knee" at ~ 4 x 10^^ eV. 

This second scenario in which extragalactic protons dominate down to the second 
knee has the following consequences: First, since the observed spectrum above the 
second knee is quite steep, oc E~^-^, the extragalactic proton flux has to cut off below 
~ 4 X 10^''eV. This can be explained as a magnetic horizon effect: Protons from 
cosmological distances cannot reach the observer any more within a Hubble time due 
to diffusion in large scale magnetic flelds jHl Ej. Second, the ankle would have to be 
interpreted as due to pair production of the extragalactic protons |H] . In particular, it 
has been pointed out recently that this model cannot afford injection of a signiflcant 
heavy component above the ankle whose secondary photo-disintegration products in the 
form of intermediate mass nuclei would produce a bump around the ankle. Injection of a 
mixed composition would also require a harder injection spectrum oc with a ~ 2.2, 
as opposed to the extragalactic ankle scenario with pure protons [8^ which requires an 
injection spectrum oc E^"^-^ . 

In the present study we point out that these conclusions can change signiflcantly 
in the presence of large scale extragalactic magnetic flelds (EGMF) which may immerse 
the UHECR sources. In fact, UHECR propagation simulations are usually performed 
for mixed compositions of nuclei without EGMF or in the presence of EGMF, but 
restricted to pure protons. This is partly due to the increased difficulty posed by nuclei 
propagating in EGMF: Especially at lower energies, deflection can be considerable and 
require considerable CPU time. As a result, this most general case starts to being 
studied only since very recently [TU[ [TT] . 

We will apply our recent work to a mixed composition similar to the one inferred for 
Galactic sources, using the EGMF obtained from the large scale structure simulations 
performed in Ref. |l2j. These EGMF are highly structured in that they reach a few 
microGauss in the most prominent structures such as galaxy clusters, but is < 10~^^ G in 
the voids. The observer is chosen in a low magnetic fleld region with some resemblance to 
Earths actual environment. We note that other simulations suggest more concentrated 
EGMF [13] and current observations do not allow to distinguish between such different 
EGMF scenarios. 

In section 2 we present our simulations, in section 3 we discuss implications for 
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the interpretation of the ankle, and we conclude in section 4. We use natural units, 
h = c = k = 1, throughout the paper. 



2. Simulations above 10^^ eV 



We restrict simulations to energies > 10^^ eV because of CPU time limitations due 
to strong deflection at lower energies and due to lack of resolution in the simulated 
EGMF. We use the same approach as described in all detail in Ref. [TU]. Specifically, we 
use a source density of 2.4 x lO^^Mpc"^, assuming the (astrophysical) sources [T^ ITH] 
follow the baryons. We inject a power law up to maximal rigidity E/Z of 4 x 10^° eV, 
whose index a can be fitted to match the observed spectrum after propagation. To take 
into account cosmic variance due to fluctuations in source positions and properties we 
simulate typically about 60 reahzations for the discrete source distribution. In addition, 
for each such realization, we dial 20 realizations for the source power which we let 
fluctuate by a factor 100 and for the injection spectral index which we let fluctuate by 
0.1 around the average fitted a, as in Ref. [Hlj. For propagation, all relevant energy 
loss processes are implemented, including pair production, pion production, and photo- 
disintegration on the combined radio, cosmic microwave, and infrared background for all 
nuclei ^Hl- In contrast to Ref. we now consider a more realistic mixed composition 
at injection that we define as follows: 
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Figure 1. Elemental abundances at a given energy per nucleon (blue asterisks) 
and resulting abundances at a given total energy for an injection spectrum oc E~'^-^ 
(histogram) . 



Let the abundance of a nucleus of atomic mass A at a given energy per nucleon be 
xa- If the spectra of all nuclei have a common slope oc the injected differential 

spectrum [duA/ dE){E) for species A at a given total energy E is then given by 

^{E)c^NxaA^-^E-\ (1) 
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where iV is a normalization constant. The modified abundance factors result from 
conversion from E /A to E. In the present work we will assume Xi = 95% for the proton 
fraction, X4 = 4% for the helium fraction, see Ref. For the remaining 1% nuclei 

with A > Awe follow Ref. [H] in adopting the distribution from Ref. JH] • The abundance 
of a stable element "^Z in the Galactic cosmic ray spectrum will be proportional to the 
product of the injection flux Eq. and the Galactic confinement time t{E/Z) which 
below the knee approximately scales as (E/Z)^^'^. Since Galactic sources are inferred 
to have an injection spectral index a ~ 2.1, Galactic cosmic ray abundances below 
the knee will be roughly oc xaA^'^ Z^'^. With the choices for xa above, this results 
in abundances at fixed total energy that are consistent with observations in the knee 
region |3JE| within a factor 2. 

Below we will find that significantly magnetized extragalactic sources should inject 
a spectrum roughly oc i?"^-^ to fit the data above 10^^ eV. The abundances xa at a given 
energy per nucleon and the ones at a given total energy resulting from Eq. (^J in this 
case are shown in Fig. ^ 

In the absence of other interactions, if nuclei of mass A at energy AE/A' 
are disintegrated into nuclei of mass A' and energy E, then we have for the fiux 
[duA' I dE) (E) of these secondary nuclei 

'"^'(E)^(i)''^(AE).N.,AA'^-^E-'. (2) 



dE ' ' \A ) dE \A' 
where in the second step we have used Eq. for the differential fiux of the parent 
nuclei. This formula holds trivially also for the primary proton flux ioi A = A = \. 




Figure 2. Realization averaged observable all-particle spectrum for injection with 
spectrum cx E~^'^ (solid horizontal line) and mixed composition as discussed in the 
text, without (red histogram) and with (blue histogram) EGMF. Cosmic variance (lines 
around blue histogram) is only shown for the case with EGMF. Also shown are some 
data from the AGASA ^5] and HiRes [201 experiments. 
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Figure 3. Average over source realizations of the all-particle spectrum (top line) and 
of spectra for different mass groups in the case of injection of a mixed composition with 
spectrum oc E^^-^ in the absence of EGMF. For particles heavier than nucleons, we 
have made the following groups for convenience: "Light elements" (blue): 2 < A < 11; 
"C group" (green): 12 < A < 24; "Si group" (black): 25 < A < 40; "Fe group" 
(yellow): 41 < A < 56. 



Photo-disintegration leads to a steepening of the spectrum compared to the 
injection spectrum. Eq. (j2)) shows that at a given energy incomplete spallation into 
nuclei A' > 1 tends to lead to secondary fluxes that are higher than for complete 
spallation into protons. The enhancement factor is given by multiplied by the 

sum of xa a over the relevant parent nuclei for A' divided by the same factor for 
protons. As a consequence, at energies around 10^^ eV where photo-disintegration starts 
to become negligible, the steepening is only significant for incomplete spallation. Fig. El 
shows that this is the case for negligible EGMF, but not in presence of considerable 
source magnetization. A similar but more moderate difference between magnetized 
and unmagnetized sources was already seen in Ref. ^0] for pure iron injection. As a 
consequence, in the absence of EGMF, a relatively hard injection spectrum oc E^"^'"^ is 
required to fit the spectrum above 10^^ eV. If we were to use a softer injection spectrum 
oc we would not fit the observed slope above 10^^ eV and we would predict a proton 

fraction around 10^^ eV of ^ 27%, too small to explain the ankle. These findings are 
consistent with Ref. P who used a homogeneous source distribution. We have verified 
that in the absence of EGMF and for the source density we use, there is no significant 
difference in the results between structured and homogeneous source distributions. 

Fig. 121 resolves the observable spectrum into individual mass group elements for an 
injection spectrum oc E'"^'"^ in the absence of EGMF. The proton fraction at 10^^ eV is 
~ 58%. 

FiglH shows that, for a fixed injection spectrum, increased spallation in EGMF can 
indeed lead to a lighter observable composition around 10^^ eV. 
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Figure 4. Realization averaged observable composition for injection with spectrum 
cx £'~^-^ and mixed composition as discussed in the text, without (red histogram) and 
with (blue histogram) EGMF. Cosmic variance (lines around blue histogram) is only 
shown for the case with EGMF. 
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Figure 5. Realization averaged observable all-particle spectrum for injection with 
spectrum oc E~^'^ (solid line) and mixed composition as discussed in the text, in 
presence of EGMF. Line key as in Fig. |21 



Enhanced, more complete spallation in the EGMF can lead to more protons which 
follow more closely the original injection spectrum because their secondary flux is smaller 
than for incomplete spallation, as explained above. This can be seen in Fig. El A similar 
tendency was already seen in Ref. JUl ioi pure iron injection. As a result, a steeper 
injection spectrum has to be chosen to fit the spectrum observed above 10^^ eV, as 
demonstrated in Fig. We note that around 10^^ eV the case of a mixed composition 
fits better than pure iron injection studied in Ref. PU] • This is because pure iron injection 
gives rise to a bump at a few 10^^ eV, somewhat flattened out by the EGMF. 
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Figure 6. Same as Fig. |S1 but for pure proton injection, with (blue histogram) and 
without (red histogram) EGMF. 
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Figure 7. Same as Fig. 13 but for injection of the mixed composition with spectrum 
oc presence of EGMF. 



The more complete spallation in case of magnetized sources consequently leads 
to a spectrum closer to the pure proton injection case, see Fig. El Note that in 
the case of pure proton injection, the observable spectrum is hardly influenced by 
EGMF. The spectrum is in fact close to the "universal spectrum" which is independent 
of the propagation mode and applies if particle number is conserved and if sources 
are distributed roughly uniformly with separations small compared the characteristic 
propagation lengths, namely energy loss and diffusion lenghts |2I]. Note that, although 
highly structured on scales < 10 Mpc, our source distribution is roughly uniform on 
larger scales, with typical source separations of ~ 30 Mpc. In contrast, the universal 
spectrum is not applicable in case of mixed compositions since the number of cosmic rays 
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is then no longer conserved. Indeed, the observable spectrum does depend considerably 
on EGMF in this case even taking into account considerable cosmic variance, as we have 
seen in Fig. |21 A numerical approach is, therefore, indispensable in this case. 



20 


1 1 1 1 Ml 1 1 1 


1 III 


1 




15 














A 

^ 10 












_ 
















5 








- 




1 1 1 




1 1 1 


1 1 1 1 



100 1000 

E [EeV] 



Figure 8. Observable composition for injection with spectrum cx E~^-^ and mixed 
composition as discussed in the text, in presence of EGMF. Line key as in Fig. ^ for 
presence of EGMF with harder injection spectrum. 

Fig. m resolves the observable spectrum into individual mass group elements for an 
injection spectrum oc E~'^'^ in the presence of EGMF. The proton fraction at 10^^ eV is 
~ 63%. 

In the presence of EGMF the composition around 10^^ eV is heavier for a steeper 
injection spectrum, see Fig. |H1 than the corresponding one shown in Fig. |3]for a harder 
injection spectrum. This is mostly due to the A""^ enhancement of heavy elements for 
a given total energy, see Eq. ([H). 

3. The Ankle 

In the previous section we saw that magnetized sources of a mixed composition require 
injection indices a ~ 2.6 to fit the spectrum above 10^^ eV, very similar to pure proton 
injection, and contrary to unmagnetized sources with a ^ 2.2. The latter spectrum 
is certainly too hard to explain the ankle and would underproduce the observed fiux 
below a few times 10^*^ eV where the injection spectrum is not modified by interactions 
anymore. In contrast, the magnetized scenario may still explain the ankle as a pair 
production feature of an extragalactic proton population dominating at that energy, as 
we are now going to argue. 

The simulations show that confinement (delay) times at energies E ~ 10^^ eV are 
~ 5 X 10^^ s, see Fig. |^ This is indeed consistent with an analytical estimate of the 
confinement time of light nuclei of energy E in large scale structures (galaxy clusters) 
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Figure 9. Histogram of delay times versus observed energy for the simulation with 
injection of a mixed composition with spectrum oc E~^'^ in the presence of EGMF. 
The contour lines are logarithmic in steps of 0.25 dex. 

of size R ~ Mpc with magnetic fields 5 ~ G, 

615(1) ~—~'"'" n^)l3M;^) 

assuming Bohm diffusion, D{E) ~ E/{3ZeB). Below 10^^ eV the simulations cannot 
be used because the Larmor radius of charged particles becomes smaller than the 
resolution of the magnetic field in the simulation. However, extrapolating Eq. Q 
to smaller energies provides us with conservatively long confinement times given that 
Bohm diffusion gives the smallest possible diffusion coefficient. Confinement times of 
light nuclei are, therefore, smaller than the age of the Universe down to ~ 10^^ eV 
and the cosmic rays of energy 10^^ eV < E < 10^^ eV will leave the strongly magnetized 
regions around the sources with roughly their original injection spectrum. Outside these 
magnetized regions space is weakly magnetized and the conditions for the universal 
spectrum should be fulfilled in this energy range [3 • In particular, the ankle should 
be reproduced as an extragalactic feature in this case, provided that the spectrum is 
sufficiently dominated by protons at these energies. 

In order to address the composition below ^ 10^^ eV, let us now assume complete 
spallation above ^ 10^^ eV, and no spallation below. The protons produced by spallation 
will pile up around 10^^ eV. Then in some range of energies around the ankle the proton 
flux will be given by the sum of Eq. Q over all A for A' = 1. On the other hand, 
the flux [duh/ dE){E) of all nuclei heavier than nucleons which are not disintegrated at 
these energies is given by the sum of Eq. over all A > 1. For the fraction of protons 
we thus obtain the rough estimate, 

{dnp/dE){E) ^ Eyi>i xa A 

{drih/dE) (E) ^ Ea>i v4°-i + Ea>i a^A ^ ' ^ ^ 
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If spallation above 10^^ eV is incomplete, the proton fraction around the ankle will 
be smaller, since the numerator in Eq. Q decreases towards Xi and the denominator 
increases due to the contribution of intermediate mass nuclei produced by spallation of 
nuclei above 10^^ eV. This tendency is confirmed by the full simulations for an injection 
spectrum oc E'"^'"^ which give a proton fraction at 10^^ eV of ~ 84% with EGMF and 
~ 58% without EGMF, compared to ~ 69% from Eq. ©. 

For a = 2.6, obtained from fitting the spectrum above 10^^ eV in the presence of 
EGMF, see previous section, we obtain ~ 44% for the proton fraction estimated from 
the simple analytical argument Eq. . The full simulations from Fig. [7| gave a proton 
fraction of ^ 63% at 10^^ eV, tending to grow with decreasing energy, and the fraction 
of the heavier mass groups is < 20%, and that of iron < 10%. Within the uncertainties 
this is probably sufficiently high to explain the ankle by pair production of dominantly 
protons. In fact, Ref. |Hj obtained that > 85% proton fraction at a given observed energy 
is needed, and from Ref. one may read off that the heavy (iron) fraction should be 
< 10%. 

Finally, at energies below the energy per nucleon where photo-disintegration lengths 
become comparable to the age of the Universe, i.e. for < 3 x 10^^ eV, the 
observed composition should approach the one injected at a given total energy. For 
the magnetized case this composition was shown in Fig. ^ as the red histogram. 

4. Conclusions 

Injection of a mixed extragalactic ultra-high energy cosmic ray population with 
composition similar to the one inferred for Galactic sources requires a relatively soft 
injection spectrum oc E^'^'^ if sources are significantly magnetized. This is in contrast 
to the case of negligible magnetic fields which requires harder injection spectra oc 
because of increased secondary fiuxes of lower energy light nuclei produced in incomplete 
photo-disintegration of heavy nuclei. Whereas the latter case requires a Galactic 
component below the ankle, the magnetized case with mixed composition may still 
explain the ankle as pair production of an extragalactic component dominated by 
protons and light nuclei piling up around 10^^ eV. This light component dominates 
the fiux around the ankle and is caused by increased spallation of the injected mixed 
composition due to considerable defiection around the sources. At energies E < 
3 X 10^^ eV the fiux may still be extragalactic, but the observed composition should 
approach the one injected at a given total energy, and may thus have a significant heavy 
component. 

If one believes that UHECR are produced by relativistic shock acceleration favoring 
spectral indices a ~ 2.2 j22], then injection of a mixed composition around considerably 
magnetized sources is disfavored. This argument may not, however, be very strong since 
there are indications that the spectral index for relativistic shock acceleration may not 
be universal 

Our results provide yet another example how structured large scale magnetic fields 
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can modify spectra and scenarios not only for discrete ultra-high energy cosmic ray 
sources [TT], but also for the cosmological sources of the diffuse flux. In addition, 
the impact of magnetic fields on the spectrum depends considerably on the injected 
composition: For pure proton injection, the effect is insignificant and the spectrum is 
close to the universal spectrum for small source separations. In contrast, for a mixed 
composition close to the one inferred for Galactic sources considered in the present work, 
magnetic fields can modify the spectrum considerably, and even stronger than for pure 
iron injection. 

We also stress that observation of a heavy composition {A) > 15 above ~ 10^^ eV 
would be almost impossible to explain by standard scenarios involving extragalactic 
sources at distances of tens of Mpc injecting compositions comparable to the ones 
inferred for Galactic sources. 
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